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Objectives 


• Look at times scales of overpressure compared to cratering. 

• What are the transient pressure and density due to atmospheric entry? 

• Does shock wave evacuate ambient gas? 

• Does transient atmospheric disturbance "settle down" during cratering? 

• Can the pressure/density field be approximated as quasi-static? 

• How does disturbance scale with impactor size? 

• What is the role of atmospheric thickness? 
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General approach 


• Perform inexpensive exploratory calculations (Phase 1). 


• Experiments to validate code and observe crater growth. 

• Follow up with more realistic coupling calculation (Phase 2). 
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MAZe Code Description 


• CFD code developed and continually evolving under various government 
contracts (DNA, NASA/JPL, NSF/USGS) 

° Formation of dust, water, ice clouds from nuclear and high explosives 

» 

° Fireball expansion, airblast propagation, boundary layer analysis 


• Relevant capabilities 

° Adaptive zoning automatically refines the grid where "needed" 

° TVD (Total Variation Diminishing) finite differencing 

Employs linearized Riemann solver, second-order accurate, flux limiting 

° Multiphase physics treats the mutual interactions between a gas mixture 
and the solid /liquid particles contained in that gas 

° ALE (Arbitrary Lagrangian Eulerian) grid techniques track 
physical interface between different type materials 
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MODELING ASSUMPTIONS 
Venus Atmosphere 

• Temperature profile from Seiff (1983) 

— 90 bar pressure, 733 K temperature at surface , 

— Treated as 100% C02 

• Non-ideal gas EOS (equation-of-state) 

t- 

— Variable y, of form P = p ( y-1) e 

— Requires y- vs-T,y-vs-e 

°111K<T< 3000 K fit from gas tables (Keenan and Kaye, 1945) 
° Assume y= 1.15 for T > 3000 K 
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Approach (Phase 1) 


• Asteroid penetration through Venusian atmosphere 

— Calculate bow shock environment ahead and behind body 

• pressures, temperatures » 

• include deceleration of asteroid 

• Treat asteroid as spherical, rigid body ^ 

— 1 km diameter, 2 gm/cc density >> Mass - 1 0 gm 

• Transparent boundary 

— Asteroid passes through planet surface without distortion 

• Minimum zone size 

— 0.125 km for 1 -km diameter impactor 
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Look at 2 impactor sizes 


• 1 km diamter and 10 km diameter. 

• Initial velocity equal to 20 km/s for each. * 

• Start leading edge at 40 km altitude . 

• 10-km body impacts at 19.5 km/s; transit time = 2.01 sec. 

• 1-km body impacts at 15.7 km/s; transit time = 2.16 sec. 

• Peak pressure is 125 kbar and 75 kbar respectively. 

• Disturbance is due to the bow wave only. 



Pressure (bar) 
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Case VEN02 1 1on diameter at 20 km/s on Venus r09 venpM 








Case VEN02 1 km diameter at 20 knVs on Venus i09venpt1 

Pressure time histories at various ranqes along around surface 
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log (R/a) log 


Vacuum or small 7T P 

Coupling Parameter: 



Crater Growth 



log (Ut/a) 


scaling laws 
C « aU* 1 



Dependence of exponents on material type 



a 

V 

Crater 

Size 

Form 

Time 

Rate of 
Growth 

Mom. Limit 

3/7 

1/3 

0.429 

0.571 

0.250 

Dry Soil 

0.507 

0.407 

0.507 

0.585 

0.289 

Wet Soil 

0.600 

0.500 

0.600 

0.600 

0.333 

Water 

0.650 

0.553 

0.650 

0.608 

0.356 

Energy Limit 

3/4 

2/3 

0.750 

0.625 

0.400 


A schematic illustration of how some of the power-law scaling relationships, which result from a 
point-source coupling parameter, are interrelated by the scaling exponent, p. 
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Vacuum crater size, time, and growth 


Impactor Size 

1-km dia 

10-km dia 

7T 2 

3.57 E-5 

157 E-4 

Final crater size 

3.9 km 

26.5 km 

R/a 

7.8 

5.3 

Time of formation 

18.2 sec 

47.5 sec 

K T = UT/a 

312 
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Growth (km, sec) 
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Pressure <KBar> Pressure <KBar> 


Test for cube-root scaling 

small 0.5% volume and large 1.5% volume 






small 1.7% volume and large 5.4% volume 
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Pressure <KBar> Pressure <KBar> 


Test for cube-root scaling (cont.) 
Need better interpolation and longer times 

small 5% volume and large 16% volume 
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small 17% volume and large 53% volume 
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Pressure <KBar> Pressure <KBar> 


Gravity-scaled crater growth 


1.2% volume 
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Gravity-scaled 
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(cont.) 
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Approach for Phase 2 
(currently underway) 

• Calculation of coupled impact, cratering, ejecta, and airblast 

— Determine partitioning of impactor kinetic energy into airblast 

— Track airblast out to times and ranges beyond crater 

• Assume impactor reaches surface intact 

• Treat body as "swarm" of various size "particles" 

— EOS includes hydrodynamic pressures under compression 

• Calculate momentum transfer to ground using ALE logic 

— Formation of impact crater and ejecta 

• EOS includes high-pressure vaporization of impactor and ground 

— Vaporized material available for airblast 
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Summary 


• This is a hard problem, all effects lead to smaller craters! 

1) Increased pressure/ drag on crater formation 

2) Slowdown of impactor due to entry drag t 

3) Breakup of impactor due to pressure distribution 

• Progress has been made in isolating variables. 

• Preliminary forms for scaling laws have be,en developed. 

• Improved calculations for transient pressure and density underway. 

• Shooting light gas gun into high pressure chamber demonstrated. 

• Equivalent explosive experiments provide supporting data. 

• Objective is to obtain useful scaling relationships for: 

Crater formation when atmospheric effects are important. 












09/199 













8182223244 PHYSICOMP 


Target EOS for Phase II 


Quartz 

• Initial Density 

• Air-filled void volume 

• Initial Bulk Modulus 

• Maximum Bulk Modulus 

• Vaporization Activation Energy 

• Initial Shear Modulus 

• von Mises Limit in Compression 

• Cohesion 

• Slope of Mohr Coulomb Surface 

• Poisson’s Ratio 


2.5 gm/cm 3 
8% 

0.45 Mb 
0.75 

1.7 * 10 10 ergs/gm 

0.173 Mb 

0.017 Mb 

50 bars 

0.75 

0.33 



